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Abstract 


To handle the complexity of modern automotive power nets, simulation-based design methods are important and suitable models of all system 
components including the battery as a main part are therefore mandatory. However, simulation models of energy storage devices are difficult to 
obtain. In particular, batteries are time-variant and strongly non-linear. An impedance-based modeling approach has been applied that copes with 
these characteristics and offers the development and parameterization of powerful models covering a wide dynamic range. As an example, this 
paper outlines the development of a NiMH battery model. Besides the impedance-based part of the model, the influences of the typical hysteresis 
effect of NiMH batteries is described in detail and an empirical modeling approach is introduced. The presented model is already successfully used 
by an automotive manufacturer which reflects the applicability of the modeling approach. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


As a successful example for impedance-based modeling and 
as an interesting application for hybrid electric vehicles, this 
paper introduces the development of a NiMH battery model. 
A non-linear impedance-based model core (see Section 2) is 
extended by a model part which considers the typical hystere- 
sis phenomenon of NiMH batteries (Section 3). This part of the 
model is parameterized in the time domain and allows together 
with the impedance-based model core for a large validity range. 
The synergies when combining impedance-based modeling 
approaches and time-domain based model parts become observ- 
able. 

The impedance-based modeling approach employs electro- 
chemical impedance spectroscopy (EIS) [1]. The measurements 
were carried out with the EISmeter, an impedance spectroscope 
that operates in galvanostatic mode and has been developed espe- 
cially for batteries and fuel cells at our institute [2]. Small-signal 
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excitations (Iac) and the evaluation of the system’s response 
enable accurate investigations at nearly any working point (state 
of charge, temperature, dc-current). The models are built by 
joining universal electrical components such as inductors, capac- 
itors and (non-linear) resistors. Physico-chemical processes can 
easily and effectively be modeled this way providing also a mini- 
mized computation effort. The method is applicable for all kinds 
of electrochemical energy storage devices [3] and [4]. 

A similar modeling approach is given in [5] for a NiMH bat- 
tery, however, without considering non-linear resistances and 
without modeling the typical hysteresis phenomenon which sig- 
nificantly influences the NiMH equilibrium-potential. Unlike 
other battery technologies, the equilibrium potential of NiMH 
batteries is not unambiguously determined by the state of charge 
(SOC) due to this phenomenon [6] and [7]. Consequently, an 
empirical hysteresis model has been developed for the pre- 
diction of the equilibrium potential and is combined with the 
impedance-based model which determines the dynamic over- 
voltages during current flows. 

The complete model implementation and a general overview 
on how to parameterize the models are given in this paper. The 
modeling approach is evaluated by comparing the measured and 
the simulated battery voltage corresponding to several current 
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Fig. 1. Impedance spectra with different superimposed direct currents [gc of a 
NiMH battery at 70% SOC and 27 °C (battery type: Sanyo Hr-DP, 6.5 Ah). 


profiles with different dynamics and current rates. All inves- 
tigations have been performed on cylindrical battery samples 
produced by Sanyo (“Hr-DP”) and Panasonic (“HHR 650 D”), 
both rated with a nominal capacity Cy of 6.5 Ah. 


2. The impedance-based core of the battery model 


Impedance spectra (Fig. 1) were obtained by applying an ac 
signal (Zac) with different defined frequencies and evaluating the 
system’s responses. After the measurement of each frequency, 
the amplitude of the ac voltage signal was automatically eval- 
uated during the measurement. For the next signal frequency 
to be investigated, the amplitude of the ac current signal was 
adapted to keep the voltage amplitude smaller than a definable 
value. Thus, the non-linearity of the battery is considered. Lim- 
ited amplitudes of approximately 3 mV per cell are usable [1]. 
For accurate investigations at different working points, bias cur- 
rents (Jac) were applied and specific cell temperatures and SOCs 
were set for the measurement of each spectrum. The temperature 
and the SOC were kept fixed as accurately as possible during the 
measurement since the impedance is a function of both states. 
To assure comparable conditions at the beginning and at the end 
of the measurement of one impedance spectrum, the maximum 
change of the state of charge during this time was limited to 
ASOC <5%. 

The following characteristics can be identified by the evalua- 
tion of the NiMH impedance spectra (Fig. 1). The spectra show 
an inductive behavior (L) at high frequencies which is mainly 
caused by the metallic connectors of the battery. The pure ohmic 
resistance (R;) can be detected by the minimum of the real part 
of the impedance spectra (at about 90 Hz). Rj reflects the limited 
conductance of the contacts, the active masses and the elec- 
trolyte. The semicircle is caused by charge-transfer processes 


Car 


Fig. 2. Equivalent electrical circuit diagram of a NiMH battery representing 
the overpotential during charging and discharging without modeling of the 
equilibrium potential. 


and charged double layers inside the battery. The changing diam- 
eter indicates the non-linearity of the charge-transfer processes. 
Different direct current rates Ia, have been superimposed to the 
small signal excitations Jac to investigate this behavior. More 
detailed information on this measurement procedure can be 
found in Refs. [3] and [4]. 

The semicircle can be modeled by a parallel connection of a 
non-linear resistor Ret and a capacitor Cg}. At lower measure- 
ment frequencies (f< 168 mHz), the diffusion behavior of the 
battery becomes apparent which is commonly modeled by a so- 
called Warburg impedance Zw. Based on this data, an equivalent 
electrical circuit of a NiMH battery can be developed (Fig. 2) 
and, additionally, be parameterized. 

Particular attention has to be paid to the non-linear behav- 
ior of the charge-transfer resistance Ret of batteries. Generally, 
this behavior can be described by the so-called Butler—Volmer 
equation with Jọ as the exchange current, n as the number of 
transferred elementary charges, 7 as the overpotential, œ as a 
symmetry coefficient, T as the absolute temperature and k as the 
Boltzmann constant (8.617 x 10E—5 eV K7!): 


I= io (ex (=) exp (=F 2) (1) 


A charge-transfer characteristic can also be observed for NIMH 
batteries. Fig. 3 exemplarily illustrates the current dependent 
values for Re, (with Ret =(dn/dD) at 27°C and 70% state of 
charge (SOC). Therefore, several impedance spectra (similar to 
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Fig. 3. Non-linear Butler-Volmer behavior (charge-transfer resistance) of a 
NiMH battery at 27 °C and 70% SOC (battery type: Sanyo Hr-DP, 6.5 Ah); black 
points: Ret values determined by a least square fit routine, black curve: pure 
Butler—-Volmer behavior (n=1, a=0.46, lọ =7 A), grey curve: Butler—Volmer 
behavior with an additional linear resistance of 0.6 mQ. 
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Fig. 4. Detailed representation of the non-linear resistor Ret (cp. Fig. 2). 


those shown in Fig. 1) have been evaluated with different super- 
imposed direct currents Jy, in a range from —80 A up to +80 A 
(~+12C rate). The evaluation of Ret has been carried out sub- 
sequently for each impedance measurement at different direct 
currents by using the equivalent circuit (Fig. 2) and a least square 
fit routine. The points correspond to the values that arise from 
this fitting procedure. A Butler—Volmer behavior is observable 
(black solid curve), however, a constant and current indepen- 
dent deviation has to be solved by adding an linear resistance. 
A microscopic physico-chemical explanation for the offset is 
unknown up to date and will be subject of further research efforts. 
However, the final result (grey curve) is absolutely satisfying 
for this kind of modeling and the later application of the NiMH 
model. 

The resistor Re can finally be modeled by means of a series 
connection of a non-linear resistor Rgy (representing Eq. (1)) 
and a linear resistor Rada as presented in Fig. 4. 

The remaining components of the electrical circuit (Fig. 2) 
have been parameterized by a least square fit routine as well. To 
enlarge the validity range to other temperatures and SOCs, the 
impedance spectra (with different Iac) have to be measured at 
various temperatures and SOCs. The presented model has been 
parameterized for an array of temperatures (— 18, 25 and 50°C) 
and different states of charge (20, 40, 70 and 90%). 


3. Hysteresis phenomenon in NiMH batteries 


For most electrochemical systems the equilibrium potential 
is unambiguously defined by the state of charge. However, the 
nickel electrode shows a significant hysteresis of the open circuit 
potential, i.e. the potential does not depend solely on the state 
of charge but also on the history of charging and discharging of 
the electrode [6,8,9]. Systematic investigations on the nature of 
the equilibrium potential hysteresis of the nickel electrode by 
Podolske Ta et al. [9] and Srinivasan et al. [8] lead to the follow- 
ing characterization of the hysteresis, that could be corroborated 
by own measurements: 


e The hysteresis is static, i.e. it remains after a charging or dis- 
charging current is switched off, even for a period of time that 
exceeds time constants of mass transport inside the electrodes 
significantly. 

e The hysteresis is independent on the preceding current rate. 

e The hysteresis is not only characterized by the boundary 
curves but also by intermediate sub-hysteresis loops. 


These characteristics of the hysteresis are in accordance with 
a domain theory, that states that the electrode material is com- 
posed of clearly separated regions (domains) each of which 
exhibits two or more metastable states ([8,9]). The actual cause 
of the metastability is yet still a matter of scientific discussion. 


OCVIV 


after continuous 


charging S 


A 


after continuous 
discharging 


0 100 
SOC I% 


Fig. 5. Schematic presentation of the hysteresis phenomenon of NiMH batteries; 
equilibrium potential is not clearly determined by SOC but rather dependent on 
pre-treatment of the battery; battery voltage is depicted without current flow 
(measurement data can be found in Fig. 11, right hand). 


One possible explanation is a change of the energy necessary for 
insertion or extraction of ions due to expansion and contraction 
of the intercalation material [9], another is a phase separation 
in the solid solution of Ni(OH)z and NiOOH in the nickel elec- 
trode [8]. Experimental results on laboratory scale electrodes 
have shown that the composition of the electrodes — especially 
the insertion of cobalt — and the ageing of the electrodes affects 
the potential level but not the shape of the hysteresis loop ([8,9]). 
However, an electrochemical explanation of the equilibrium 
potential hysteresis is out of the scope of this paper; in the 
following a model will be developed that describes the effect 
quantitatively, can be parameterized by electrical measurements 
and is efficient with respect to computation time in a computer 
simulation model. 

The lower and an upper-boundary equilibrium-potential 
curve measurable at open circuit conditions after charging or dis- 
charging periods is shown in Fig. 11 (measurements) and Fig. 5 
(schematic illustration). Moreover, the intermediate range can be 
reached if the battery is operated with shallow cycles (depicted 
with small cycles). If the battery is for example continuously dis- 
charged until “A” and then re-charged, the equilibrium potential 
(measured without current flow) will follow the upper part of 
the small cycle (“A to “B”) and then finally describe the outer 
curve until 100% SOC. The necessity of an accurate modeling 
of this effect becomes obvious in Fig. 6 which illustrates that the 
hysteresis effect can be even more significant than the dynamic 
overvoltages up to absolute current rates of |Jpatt| < 20 A. 

According to systematic investigations, the hysteresis effect 
depends exclusively on the charged or discharged amount of 
ampere hours. Comprehensive tests and conclusions showed 
that the current rate does not influence the equilibrium-potential 
level. A test with alternating discharging and charging pulses 
(different current rates but constant amounts of 3000 As) was 
carried out (Fig. 6). Rest periods of 1000s duration were per- 
formed after each pulse to check the equilibrium potential. The 
gap between the straight lines (upper and lower equilibrium- 
potential curves) is not affected by the current rate which is in 
line with the findings in [6]. 
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Fig. 6. Alternating charging and discharging pulses with different current rates 
but with a constant amount of 3000 As per pulse and rest periods of 1000s in 
between; the hysteresis effect (gap between the straight lines) is independent of 
the current rate (battery type: Sanyo Hr-DP, 6.5 Ah/70% SOC, 27°C). 


3.1. Modeling of the hysteresis phenomenon 


Generally, the equilibrium potential (Up) of a battery cell 
can be modeled by a SOC-dependent voltage source. Since the 
NiMH equilibrium-potential behavior is seriously influenced by 
the above described phenomenon, the voltage source has to be 
adapted to a more complex model. The complete equivalent- 
circuit is depicted in Fig. 7, the detailed algorithm to calculate 
the equilibrium potential Uo(SOC, hyst) is addressed below. 

Uo is calculated by the sum of the lower SOC-dependent 
boundary voltage ULg and an additional voltage Unys, which 
considers the hysteresis effect of partial cycling (Eq. (2)). The 
maximum value of Uo is limited by the upper boundary Uyp, 
the voltage Unyst is limited SOC-dependently by the gap of the 
lower and upper boundary 


with Up(SOC) < Uyp(SOC) 
(2) 


Both, U_p(SOC) and Uyp(SOC) can easily be modeled by a 
fixed polynomial approximation (5th order) with the SOC as the 
input value. The SOC is calculated according to Eq. (3) with 7 
as the battery current and Cy as the nominal battery capacity 
(convention: charging J > O/discharging 7< 0). The gassing reac- 
tion can be taken into account by including the gassing current 
Igas When integrating the battery current as shown in Eq. (3). 
The gassing current is a non-linear function of temperature and 
potential. A proposal for an appropriate modeling approach of 


Uo(SOC) = ULB (SOC) + Uhyst 


U,(SOC, hyst) 


Cy 


Fig. 7. Equivalent electrical circuit diagram of a NiMH battery considering the 
equilibrium potential which shows a hysteresis effect. 
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Fig. 8. Closer examination of the shape of the inner course of the equilibrium 
potential between the outer boundaries of the hysteresis curves (cp. Fig. 5). 


the gassing reaction is given below and its necessity is discussed. 


t 
Sizo = Ieas)dt | o 


SOC = SOC(t = 0) + 
CN 


with SOC < 100% (3) 


Unys has to be calculated during the simulation by a more com- 
plex algorithm. Therefore, a new parameter Qhyst has to be 
defined (Eq. (4)) 


Gi J Tir with O< Ora < Qnyst_max (4) 


Qhyst Sums up the charge-transfer current during a microcy- 
cle of the hysteresis effect. Qpyst is equivalent to 0 when the 
battery only has been continuously discharged and the equi- 
librium potential follows the lower boundary (see Fig. 5). 
In case the battery is subjected to a microcycle (e.g. from 
“A” to “B”, Fig. 5), Qnyst increases and is finally limited at 
the value Qnyst_max when the upper boundary is reached (at 
“B”). Qhyst max has to be determined by measurements (see 
Section 3.2). Unyst is consequently equal to OV at Qnyst=0 
As and Uhyst = Uug(SOC)—-ULB(SOC) = Uhyst_max(SOC) at 
Qhyst = Qhyst_max- The maximum hysteresis voltage Ubyst_max has 
also to be deduced from measurements (see Section 3.2). 
Considering these basic principles, only the shape of the 
course of the voltage during a shallow cycle between the outer 
boundaries has to be considered (Fig. 8). Therefore, the maxi- 


7 : is h 
mum stretching of the inner voltage characteristic DU oct max 


for charging and an 3 for discharging (at Qnyst_max/2) 
have to be evaluated by measurements (see Section 3.2). Since 
the hysteresis voltages at Qhyst=0, at Qnyst = Qhystmax/2 and 
Qhyst = Qhyst_max are thus determined, both curves can be approx- 
imated by polynomial equations of second order. By means of 
the limitations Unyst_max(SOC) and U_p(SOC), the stability of 
the hysteresis model over the complete SOC range is assured. 
Completely reasonable voltage values cannot be achieved since 
the lower and the upper boundaries cannot be overstepped. 
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Fig. 9. Measurement of the inner voltage curves between lower and upper 
boundary, two different discharging periods (B: after 0.7 Ah charging, C: after 
1.2 Ah charging) (battery type: Sanyo Hr-DP, 6.5 Ah/70% SOC, 27°C). 


A more dynamic operation with smaller Ah throughputs than 
Qhyst.max leads to processes inside the two already illustrated 
(Fig. 8) voltage curves. The measured data in Fig. 9 illustrate an 
appropriate test: 

The battery has been continuously discharged to 70% SOC 
(which means Qhyst equals 0 Ah). Subsequently, the battery has 
been charged with 10 A up to Qnyst = 0.7 Ah (9.3% Cn, point “B” 
in the graph) and charged to Qhys = 1.2 Ah (18.5% Cy, point 
“C”) respectively. Rest periods (1000s) have been introduced 
to determine the equilibrium potential and Unys, during the test. 
After reaching point B and C respectively, the battery has been 
discharged with —10 A until Qnyst =0 Ah again. A very similar 
symmetry of both curves (“A —> B —> A” and “A —> C > A”) can 
be observed considering a scaling factor which is linear to the 
actual amount of Qhyst. This correlation can be used for modeling 
(Eqs. (5) and (6)) 


discharge _ discharge Qhyst 
A Unyst =A Unyst_ max Q (5) 
hyst- max 
charge _ charge Ohyst_ max — Ohyst 
A Unyst =A U yst- max (6) 
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Fig. 10. Exemplary illustration of the modeling of Unsyt(Qnyst) when the battery 
charge throughput is smaller than Qnyst_max- 


The outer voltage curves (Fig. 8) can be approximated by sec- 
ond order polynomials as mentioned. Using Eqs. (5) and (6), 
the same procedure can be applied when the Ah throughput is 
smaller than Qhyst max. To explain the employed algorithm, a 
charging/discharging sequence is explained by means of Fig. 10 
(schematic, no measurement) as an example. 

The sequence starts at point “A”. Initially, the sign of the cur- 
rent is determined. Since the first pulse is a charging pulse, the 
voltage starts to rise according to the charging curve “‘ACB”. 
When charging is finished at some point (here: “C”’), this point 
corresponds to the first point of the new discharging curve 
“CDA”. The second point is the starting point “A”. The third 
point (“D”) can be determined by considering half of the charged 
value (charging from “A” to “C”) and using Eq. (5). This is suf- 
ficient for the calculation of the polynomial coefficients for the 
discharging curve “CDA”. 

Instead of a finished discharging to point “A”, the discharging 
may be finished in point “E” and charging starts again. The first 
point for the calculation of the next charging period is point “E”, 
the second is again the saturation point B. For the calculation of 
the charging curves, point “B” is always relevant as the reference 
point, and for the calculation of the discharging curve point “A” 
is the reference point. The still needed third point for the coef- 
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Fig. 11. (left hand) Measured equilibrium-potential values of a NiMH battery (points); lower curve measured after discharging steps (—10 A) with rest periods 
(1000 s) in between, upper curve after charging steps (+10 A); (right hand) Measurement data of the inner voltage curves according to the measurement procedure 


illustrated in Fig. 9 (battery type: Panasonic HHR 650D, 6.5 Ah/27 °C). 
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ficients of the new charging curve (“EFB”) can be calculated 
again by using the similarity (Eq. (6)). 


3.2. Parameterization of the hysteresis phenomenon 


In the following paragraphs details are given on how to 
parameterize the model for the hysteresis. 


3.2.1. Determination of Utp(SOC) and Uyp(SOC) 

Measurement results to parameterize the lower and the upper 
boundary are illustrated in Fig. 11 (left hand). Therefore, a fully 
charged NiMH battery (nominal capacity 6.5 Ah) has been dis- 
charged stepwise with a current rate of 7=—10A. Again, rest 
periods of 1000 s after each discharging event (10% of Cyn) have 
been made to measure the equilibrium-potential. Subsequently, 
the battery has been fully charged (J= 10 A) using the same pro- 
cedure. For the SOC determination, coulombic losses due to 
gassing have been considered by means of the proposed gassing 
model (Section 4). For the implementation, the boundary curves 
are approximated by polynomials of fifth order. 


x š h 
3.2.2. Determination of Qhyst max, AU 8. and 


hyst- max 
discharge 
A U hyst- max 


The procedure according to Fig. 9 has been applied at three 
different SOC values (initiated at 35, 60, and 80% SOC). 
Measurement data can be found in Fig. 11 (right hand). The 
deviations between the lowest points of the inner curves and 
the lower boundary curve (ULp(SOC)) are addressed in the 
discussion part. Saturation values of Qhyst.max % 0.8. ..1.2 Ah 
have been found. The maximum stretching of the inner voltage 
curves has been IE ine ate x~ 17...20mV and a x 
5...10mV. The stretching during charging was always signifi- 
cantly higher than in case of discharging. The values have been 
determined for cylindrical NiMH batteries with Cy =6.5 Ah 
(samples from Sanyo and Panasonic). 


3.2.3. Determination of Unyst max 

The maximum hysteresis voltage Uhystmax(SOC) corre- 
sponds to the difference of Uyg(SOC) and U_p(SOC). Values 
of up to 70 mV have been found 


Uhyst-max (SOC) = Uup(SOC) — Urg (SOC). (7) 
4. Proposal of a gassing model 


Measurements were performed at different battery temper- 
atures to evaluate the gassing characteristic. Therefore, a fully 
charged cell was overcharged with different constant currents in 
arange of 0.325 A (0.1 C) up to 3.25 A (1 C). The corresponding 
cell voltage was collected after its stabilization. Fig. 12 gives the 
measurement data by using a semi-logarithmic illustration. 

The lines represent the results of Eq. (8) which can be used as 
an additional model part. T corresponds to the battery temper- 
ature and U to the battery voltage, however without the ohmic 
overpotential which does not affect electrolyte decomposition. 
Uo corresponds to the decomposition voltage of water with 
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Fig. 12. Gassing currents at 30°C, 2°C and —28 °C and different cell voltages; 
the temperature of the battery could be regulated to a range of +2 °C around the 
illustrated values; the lines represent the results of the proposed gassing model. 


1.23 V. The remaining parameters are freely chosen values; the 
values listed in Table 1 result an accurate match of the measured 
data (see Fig. 12) 


Igas = Ipexp(a(U — Up) )exp(A(T — To)) (8) 


The necessity of this modeling part depends on the application 
of the model. If the voltage behavior has to be investigated for 
quite short periods, the gassing model can be neglected probably. 
This is exemplarily demonstrated by Fig. 13. The left hand figure 
illustrates the simulated voltage responses to a profile with alter- 
nating current steps (+30 A, ~4.6C) at 27°C and 70% SOC. 
Rest periods are inserted between the charging and discharging 
steps. The duration of each step and the rest periods amount 
to 30s. The simulations have been performed with and without 
gassing model. Even after a period of half an hour, no voltage 
deviations can be observed between both approaches. 

However, the gassing influences cannot be neglected if the 
simulation of the battery efficiency is aimed. This is illustrated 
exemplarily by means of the right hand figure. The charging cur- 
rent is reduced by the coulombic losses (gassing current) which 
lead to a clearly decreasing state of charge (~ 1.5% SOC during 
this test). This effect intensifies with more intensive charging 
cycles, with increased SOC and with increased temperature. 
However, since the description of the pure voltage behavior is 
the main focus of this paper, the gassing influences were not 
considered for the simulated (quite short) voltage courses given 
in Section 5 (Igas = 0). 


Table 1 
Parameters for modeling of the gassing behavior by means of Eq. (8) corre- 
sponding to measurement data illustrated in Fig. 12 


Ip 1.808 mA 
To 25°C 

a 25.74 V7! 
B 0.078°C7! 
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Fig. 13. Simulation with and without gassing model of alternating current steps at 27 °C with 30s duration, +30 A battery current and rest periods in between; the 
gassing influences on the voltage behavior are almost negligible for the simulated period of 0.5 h (left hand), the influences on the SOC calculation are not negligible 


(right hand). 
5. Model evaluation 


Comparisons among measured and simulated voltage curves 
over time are presented in this section. Step profiles with differ- 
ent dynamics (respectively step duration) and a wide range of 
current rates were applied to a battery (Sanyo Hr-DP, 6.5 Ah). 
Fig. 14 exemplarily illustrates the procedure of the current steps 
(left hand) and the corresponding cell voltage (right hand). It has 
to be noted that the cell voltage changes even though no current 
is applied. The recovering voltage is caused by the double layer 
capacitance Cq and the diffusion processes of the charge car- 
riers (see Section 2). Figs. 15—17 illustrate the simulation and 
measurement results according to this kind of input profile. 

Fig. 15 gives the results for the steps with 10s duration. The 
test was carried out at 27°C and 20% SOC. The current range 
amounted to 0.5 to 15 times the C-rate. 

A further illustration compares measured and simulated data 
by using shorter steps of 1s duration each. This test was per- 
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formed at 50°C and 20% SOC. The current range corresponds 
to the current range which was the basis for the data in Fig. 15 
as well. 

Fig. 17 depicts a comparison for input currents in a range 
from 0.5 to 4.6 times the C-rate at 27°C and 70% SOC. Alter- 
nating current pulses with a lower dynamic of 100s duration 
were applied. 

Finally, a current profile is illustrated that forces the hystere- 
sis effect and illustrates the advantages of the actual hysteresis 
model. Therefore, current pulses with 10 A were applied, how- 
ever with different duration and therefore different amounts of 
charge. The current profile is given in (Fig. 18). The cumulative 
amount of charge of the three consecutively repeated discharge 
pulses corresponds to the amount of Ah of the preceding charg- 
ing step. Moreover, the total amount of charged Ah was increased 
subsequently. Hence, hysteresis curves of different sizes (cp. 
Fig. 9) have to be simulated accurately to match the measured 
voltage (Fig. 19). 
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Fig. 14. Exemplary illustration of input current steps used for model evaluation; repetitions of these current steps have been applied to a real battery and to the model 


with different step durations; results are depicted subsequently. 
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Fig. 15. Comparison of measured and simulated battery voltage corresponding 
to an alternating profile with absolute battery-currents of up to 100 A at 20% 
SOC and 27 °C (battery type: Sanyo Hr-DP, 6.5 Ah). 
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Fig. 16. Comparison of measured and simulated battery voltage corresponding 
to an alternating step profile with absolute battery currents of up to 100 A at 20% 
SOC and 50°C (battery type: Sanyo Hr-DP, 6.5 Ah). 


6. Discussion 


An accurate agreement between the simulated and the mea- 
sured battery behavior has been demonstrated. These results 
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Fig. 17. Comparison of measured and simulated battery voltage corresponding 
to an alternating profile with absolute battery currents up to 30 A at 70% SOC 
and 27 °C (battery type: Sanyo Hr-DP, 6.5 Ah). 
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Fig. 18. Input current profile to force the hysteresis effect, the discharged amount 
of each discharge triple corresponds to the preceding charging step. 


have been achieved by combining an impedance-based model 
core which describes the battery dynamic in a wide range and 
a hysteresis model that successfully describes the hysteresis 
of the equilibrium-potential of NiMH batteries. However, the 
hysteresis model is a purely phenomenological approach and a 
certain error between the modeling assumption and the mea- 
surement data becomes obvious. For modeling, small inner 
curves with connecting points to the outer hysteresis curve have 
been assumed (cp. Fig. 5 (schematic)). This approach has been 
accepted even though measurements indicate a voltage error of 
approximately 10 mV where the inner curves do not really join 
the lower boundary curve (cp. Fig. 11 (measurement)). More 
theoretical explanation for the hysteresis phenomenon would be 
beneficial for the understanding and model implementation. The 
hysteresis phenomenon of NiMH batteries has been previously 
investigated and explained on a qualitative basis (e.g. in [6,8,9]), 
but no consistent model has been presented that explains the 
effect quantitatively. 
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Fig. 19. Voltage response to the current profile which is depicted in Fig. 18 at 
27°C and 70% SOC (battery type: Sanyo Hr-DP, 6.5 Ah). 
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Measurements have been shown for 27°C and 50°C within 
this paper. The impedance based modeling approach has been 
proven down to —18°C as well. However, the collection of 
impedance data at deep temperatures is more difficult (the 
model parameterization, respectively). On the one hand, the 
battery heats up quite fast during tests with bias current. On 
the other hand, the transport of charge carriers becomes more 
relevant with decreased battery temperature. Investigations on 
transport processes however demand quite high Ah throughputs 
which violate the rules of impedance spectroscopy. “Quasi- 
steady” conditions with no significant changes of the SOC have 
to be kept during the measurement. Consequently, the purely 
impedance modeling approach is restricted. An enlargement 
of the model’s validity range (including large Ah-throughputs 
combined with deep temperatures) could be achieved by addi- 
tional model parts such as a separated transport model which 
demands no impedance measurements. This has been applied 
for the modeling of lead-acid batteries already [10]. 

Furthermore, the battery model has to be equipped by an 
appropriate gassing model if the simulation of the battery effi- 
ciency is focused (not only the voltage behavior for a short 
duration). A possible gassing model has been proposed which 
can be connected in parallel to the RC-element, the Warburg 
impedance and the equilibrium-potential calculation (Fig. 7). 
This approach has been employed successfully as well for lead- 
acid battery models [11]. 


7. Summary 


An impedance-based modeling approach, that copes with the 
specific battery characteristics and offers the development and 
parameterization of powerful models covering a wide dynamic 
range, has been applied to NiMH battery technology. This 
modeling approach has been successfully applied to further bat- 


tery technologies such as Li-ion, lead-acid and electrochemical 
double-layer capacitors [4]. Besides the impedance-based core 
of the model, the typical hysteresis in the equilibrium voltage, 
which is very pronounced for NiMH batteries has been described 
in detail and an appropriate modeling approach has been intro- 
duced. Finally, the NiMH battery model has been evaluated with 
comparisons among measured and simulated voltage data in the 
time domain by using input current profiles at different states 
of charge, with significantly different current rates and vari- 
ous time constants (Section 5). Limits and valuable upgrades 
of the NiMH battery model have been finally discussed. The 
proposed gassing model has to be considered if the simulation 
of the battery efficiency is focused. 
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